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We introduce a simple prescription for calculating the spectra of thermal fluctuations of 
temperature-dependent quantities of the form STit) — J d^fST{f,t)q{r). Here T{f,t) is the lo- 
cal temperature at location r and time t, and q{f) is an arbitrary function. As an example of a 
possible application, we compute the spectrum of thermo-refractive coating noise in LIGO, and find 
a complete agreement with the previous calculation of Braginsky, Gorodetsky and Vyatchanin. Our 
method has computational advantage, especially for non-regular or non-symmetric geometries, and 
for the cases where g(r) is non- negligible in a significant fraction of the total volume. 



I. INTRODUCTION AND MAIN RESULTS 

The theory of time-dependent thermodynamical fluc- 
tuations has been extensively developed for the past cen- 
tury. One of the fundamental results in this field is 
the Fluctuation-dissipation Theorem (FdT), which was 
originally formulated by Callen and Welton in 1951 [l|. 
Several different formulations of the FdT have been in- 
troduced since then, see [2,] for a review. In this short 
paper we develop a formulation which is suitable for cal- 
culating the spectra Sgj,{uj) of thermal fluctuations of 
temperature-dependent quantities of the form 



(1) 



ST = / (Fr5T{r,t)q{r). 



Here 5T{r, t) is the local temperature at location r and 
time t, and q{f) is an arbitrary function. 

Computations of this kind may be relevant for the de- 
sign of interferometric gravitational-wave detectors like 
the Laser Interferometer Gravitational Wave Observa- 
tory (LIGO). Random thermal fluctuations are expected 
to be the dominant noise source for LIGO at frequencies 
between 10 and 100 Hz; see, e.g., [s*]. This noise is mostly 
due to the thermal motion of the LIGO mirror surfaces, 
and has been extensively studied both theoretically and 
experimentally by several groups around the world; sec, 
e.g., [3l and references therein. More recently Braginsky, 
Gorodetsky, and Vyatchanin (Q, hereafter BGV) have 
identified a different kind of thermal noise, which they 
called the thermo-refractive (TR) noise. The TR noise 
can be understood as follows: 

Consider a laser beam which passes through (the part 
of) one of LIGO test masses. The index of refrac- 
tion of the test-mass material is strongly temperature- 
dependent. Thus local thermodynamical fluctuation in 
the temperature ST{r) result in the fluctuation of the 
overall phase of the laser beam, as measured by an inter- 
ferometric experiment. For small expected temperature 
fluctuations, the variation of the phase is of the same 
form as in Eq. ([T]). Thus a computation of the thermo- 
refractive noise amounts to that of the fluctuations in ST. 
BGV have computed the TR noise for a plane-parallel 
geometry of the thin mirror coating. Their computation 



was based on the multi-dimensional Langevin equations 
and involved a calculation of the spacial temperature cor- 
relation functions. While BGV approach works fine for 
simple geometries, it becomes numerically cumbersome 
in more general geometries, i.e. in a case when the light 
goes through one of the test masses and the light-beam 
radius is comparable to that of the test mass. 

By contrast, our computational approach will be based 
on a direct application of the FdT, in the spirit of our 
earlier treatment of the mechanical thermal noise ([^, 
hereafter L98). We will show that in order to compute 
the fluctuations in T, one needs to perform the following 
mental experiment consisting of 3 steps: 

1. Periodically inject entropy into the medium, with 
the volume density of the entropy injection given by 



Ss(f) 
dV 



FQCOs{ujt)q{f), 



(2) 



where Fq is an arbitrarily small constant. 

2. Track all thermal relaxation processes in the system 
(e.g., the heat exchange between different parts of the 
system) which occur as a result of the periodic entropy 
injection. Calculate the total entropy production rate 
and hence the total dissipated power Wdiss which occurs 
as a result of the thermal relaxation. 

3. Evaluate the spectral density of fluctuations in ST 
via the formula below: 



(3) 



cf. Eq. (1) of L98. Here Lj^2nf. 

In the following section we give a proof of the above 
prescription. In section 3, as an illustration of the 
method, we compute the thermo-refractive coating noise 
in the BGV geometry. We conclude with some brief gen- 
eral discussion in section 4. 



II. PROOF 

It may well be possible to verify Eq. ^ by appeal- 
ing to one of the existing formulations of the FdT, which 
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treat the fluctuations of generalized thermodynamic vari- 
able (see [2]). However, the gravitational- wave commu- 
nity (including this author) is much more familiar with 
the FdT for a generalized mechanical coordinate of the 
system, as given by Callen and Welton in 1951 fl!]. We 
thus think it is instructive to construct a proof of Eq. ^ 
using a mechanical coordinate. 

For this, we mentally introduce a set of non-intrusive 
mechanical thermometers into the system. Our ther- 
mometers are localized ensembles of identical harmonic 
oscillators, which are assumed to 

(a) be sufficiently densely packed into the system (we will 
make this more precise shortly), 

(b) have total mass and heat capacity which are vanish- 
ingly small compared with those of the original system 
(for the latter it is sufficient to assume that the number 
of thermometers is much smaller than the number of par- 
ticles in the system), 

(c) have a proper angular frequency loq which is much 
higher than the angular frequencies 27r/ at which Sgrp{f) 
is computed, and 

(d) be thermally coupled to the system on a timescale 
much shorter than 1//. 

Consider now an operator 



(4) 



where Xi is the displacement of the i'th oscillator. For 
densely packed oscillators, the equation above can be 
written as 



Ti = / (frn{f) < x2(f) >, 



(5) 



where n{f) is the spacial density of the oscillators, and 
< a:^(r) > is the average x"^ at a radius r. Note that as 
the number of oscillators increases, the Eq. ([5]) becomes 
better defined and more precise. For sufficiently dense 
packing, 



< x'^{r) >: 



rT(r), 



(6) 



where fcs is the Boltzmann constant, and m is the oscil- 
lator mass. Here, when we say "sufficiently dense", we 
mean that in a minimum volume for which the notion 
of local temperature T(f) is meaningful, there should 
be a number of oscillators ^1. If such volume has 
Ni ^ 1 particles, then it should have N2 oscillators 
where TVi > iV2 > 1. 

We now choose the density of the oscillators be 



n(r) = q{r), 

Kb 



(7) 



where q(r) is the form factor from Eq. ([T]); we can always 
rescale q{r) so that the oscillators are sufficiently densely 
packed. With this choice of n(r), the dynamical variable 
Ti closely tracks the thermodynamical variable T. 

We can now use the Callen- Welton formulation of 
the FdT to find the fluctuation Sf{f) in % We follow 



closely the steps described in L98: 

Step 1. We introduce a periodic perturbation of the 
form of the interaction Hamiltonian 



i?i„t = --FoCOs(27r/i)Ti, 



(8) 



and consider response of the system to this perturbation. 
From Eq. ^ we see that physically such perturbation 
amounts to a periodic change in the rigidity of each of 
the oscillator, or, equivalently, in a periodic change 5ujq 
in the oscillator proper frequency: 



5wo 



Fq cos(27r/t) 



(9) 



If the oscillators were not thermally coupled to the sys- 
tem, their energy would track adiabatically the change in 
the proper frequency, SE/E = Sujo/loq. Once the thermal 
coupling is included, their energy change is given by 



SE^ 
E, 



E,, 



(10) 



where dQi is the energy input from the I'th oscillator to 
the local thermal bath. By construction the oscillators' 
thermal coupling to the system is much more rapid than 
1//, and thus on average their energy does not change (it 
remains ksT). Therefore, once averaged over the local 
volume dV, Eq. pil)l gives 

5Q = ^££^£(|LZ!) X ksTn{r) x dV. (11) 



We now substitute Eq. ([7]) into the above equation, and 
get 



5s 
dV 



TdV 



Fo cos{2T:f t)q{r). 



(12) 



This is the density of the local entropy injection as a 
result of the periodic perturbation driven by Hi^f 

Step 2. We compute the total power Wdiss which 
is dissipated in the system as a result of the periodic 
forcing. This means that we compute the total entropy 
production in the system; from Eq. (jl2[) we see that it is 
zero to first order in Fq. Not so to the second order: the 
entropy injection leads to temperature inhomogeneities 
which in turn lead to thermal relaxation processes, i.e. 
to the heat exchange between different parts of the 
system. The thermal relaxation leads to the net entropy 
production rate which is second order in Fq. 

Step 3. We compute Sgj.{f) by using Eq. This 
completes justification for and physical description of the 
computational procedure outlined in the introduction. In 
next section we work out a practical example, as a useful 
crosscheck of our approach. 
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III. THERMO-REFRACTIVE NOISE IN LIGO 
MIRROR COATING 

BGV studied the thermo-refractive noise in the thin 
optical coating of LIGO mirrors. We re-derive their result 
using the direct approach developed above. 

As BGV have explained, the relevant variable for the 
coating thermo-refractive fluctuation is 



6T 



1 



'(x^W}/rlp-z/l 



dxdy / dz5T{r,t)e' '-- "'"e 

0' J -aa Jq 



(13) 

where tq is the effective beam size, / is the effective coat- 
ing thickness, x and y are the rectangular coordinates 
along the mirror face, and z is the coordinate along the 
beam, chosen so that z = at the coating outside bound- 
ary. Within our approach, we need to, as a thought ex- 
periment, inject an oscillating entropy perturbation of 
the form in Eq. fj]), where 



q{f) 



1 



27^ 



(14) 



The heat injection leads to periodic temperature pertur- 
bation ST{'r), which we compute below. We can then 
find the dissipated power [see Eq. (35.1) of Landau and 
Lifshitz [Z] and Eq. (5) of 



J {{V5Tf) , 



(15) 



where k is the thermal conductivity, and < ... > stands 
for averaging over the oscillation period 2-k /lj. 

As B GV have noted, the thermal diffusion lengthscale 
Itii = a/ k/ {Cpuj) satisfies the following inequality: 



ro > hh > I- 



(16) 



Therefore, (1) the heat diffusion is almost exclusively in 
the z direction, and (2) we can consider all of the entropy 
to be injected at the outer coating boundary, with the 
surface density 



dA 



(17) 



The temperature perturbation satisfies the diffusion 
equation 



d6T _ d^ST 
and the boundary condition 



fddT\ ^ d_5s_ 



(18) 



(19) 



The latter simply states that all the heat injected at the 
boundary is transported inwar ds[g. Here D = n/{pC), 
where p is the density and C is the heat capacity. 
It is now straightforward to find 5T and its gradient: 



dT TFqlu 




z sm 




z — Lot \ e 



dz KTxrk \ y 2D ' J V V 21? 
We can now find the dissipated power: 

F^Lv^T ID 



(20) 



,1^ (dST\ 
T\\dz , 

Finally, from Eq. we get 



47rrQK V 2aj 



(21) 



SsAf) = 



Trr^y/LuCpn 



(22) 



This resuh is identical to Eq. (9) of BGV. 



IV. DISCUSSION 

Our direct approach to the thermo-refractive noise 
has 2 advantages, as compared to the previously used 
methods: 

1. It is computationally easier, especially when the 
form-factor q{r) is significant over a large fraction of the 
system's volume. For example, a broad beam of light 
passing through a LIGO test mass will be sensitive to 
the thermo-refractive fiuctuations in a large fraction of 
the test mass' volume. 

2. In elucidates geometries and material properties 
which may make the thermo-refractive noise high. The 
thermal relaxation which is induced in the process of our 
mental experiment depends on the spacial variation of 
the form- factor q{'r) and (less obviously) on that of the 
heat capacity. 

Future work may include unified simultaneous modelling 
of thermo-refractive and thermoelastic noises, since they 
both originate from T-fluctuations and may be highly 
correlated. 
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